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ABSTRACT 

We present a catalog of radio afterglow observations of gamma-ray bursts (GRBs) over a 14 year 
period from 1997 to 2011. Our sample of 304 afterglows consists of 2995 flux density measurements 
(including upper limits) at frequencies between 0.6 GHz and 660 GHz, with the majority of data taken 
at 8.5 GHz frequency band (1539 measurements). We use this datasct to carry out a statistical analysis 
of the radio-selected sample. The detection rate of radio afterglows has stayed unchanged almost at 
31% before and after the launch of the Swift satellite. The canonical long-duration GRB radio light 
curve at 8.5 GHz peaks at 3-6 days in the source rest frame, with a median peak luminosity of 10"^^ 
erg s~^ Hz""'^. The peak radio luminosities for short-hard bursts. X-ray flashes and the supcrnova- 
GRB classes are an order of magnitude or more fainter than this value. There are clear relationships 
between the detectability of a radio afterglow and the fluence or energy of a GRB, and the X-ray 
or optical brightness of the afterglow. However, we find few significant correlations between these 
same GRB and afterglow properties and the peak radio flux density. We also produce synthetic light 
curves at centimeter and millimeter bands using a range of blastwave and microphysics parameters 
derived from multiwavelengtli afterglow modeling, and we use them to compare to the radio sample. 
Finding agreement, we extrapolate this behavior to predict the centimeter and millimeter behavior 
of GRBs observed by the Expanded Very Large Array (EVLA) and the Atacama Large Millimeter 
Array (ALMA). 

Subject headings: gamma rays: bursts — cosmology: observations — hydrodynamics — radio continuum; 
general 



1. INTRODUCTION 

Our understanding of gamma-ray bursts (GRBs) has 
advanced rapidly since the discovery of the long-lived 
"afterglows" at X - ray, optical and radio wavelengths 
(iCosta et all 119971: Ivan Paradiis et all Il997t iFrail et al.l 
|1997[) . Since that time many more 7— ray and hard X- 
ray satellites have monitored the sky for GRBs providing 
accurate localizations and enabling detailed followed up 
by the space and ground-based facilities at longer wave- 
lengths. As a result of this large and sustained effort, 
several hundred afterglows have been detected over the 
last 15 years. The afterglow study has exposed a rich di- 
versity of GRB behaviors a nd several distinct progenitor 
classes (|Gehrels et al.|[2009t) . While many of the early ad- 
vances in the field were the result of studying individual 
GRBs and their afterglows, the field has now matured 
to the point where much more can be learned from the 
study of large samples. 

Statistical analysis have long been carried out 
on GRB-only sampl es (e.g. iKouveliotou et all 119931 : 
ISakamoto et al.ll201lD but including afterglow data has 
only been done recently. Compilations of X-ray and 
optical light curves from the Swift satellite have re- 
vealed complex bu t canonical light curve behaviors 
(iNousek et al.l 120061: iZhang et al.l l2006t [Melandr i et al.l 



terglow luminosities (iLiang fc Zhan^ l2006t iKann et "all 
l200alNardini et al.ll2006rr 

Comparative studies of GRBs and their afterglows 
have proven even more useful. The standard fireball 
model has been tested through spectral and tcmpo- 
ral comparisons of X-ray and optical light curves (e.g. 
lOates et al.ll2009l l20nl: ISchulze et al.ll201lD . Significant 
differences have been found between the mean brightness 
and the reds hift distrib ution of the Swift and the pre- 
Swift GRBs (iBerger et a l. 2005h [Jakobsson et al.l [20061: 
IKann et al.ll2010D . Correlations have been found for both 
short and long-duration bursts between the gamma-ray 
fluence and the X-ray and optical afterglow brightness 
(iGehrels et al.ll2008HNvsewander eraIll2009l:lKann et al.l 
I2OIII). A populatio n of dark bursts has been identi- 
fied (iRol et aLll2Q05|) . as has a population of nearby, low 
luminosity events ( Sazonov et al.ll2004l: [Soderberg et al. 
2Q06al: IPian et al.ll20Q6l: ILiang et al.l 120071: IKann et al 
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120081: IRoming et all [20091: lEvans et'aLl [20091). Optical 
catalogs have d erived mean dust extinction laws of GRB 
host galaxies (jKann et al.l 120061: [Starling et all I2007D . 
and have lead to claims of clustering of the optical af- 
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In comparison, very little effort has gone into compiling 
radio afterglow data and carrying out correlative studies 
([Berger 2004a). A catalog of the first five years of radio 
afterglow data was produced by iFrail et all (|2003l ) , but 
comparisons were limited just between detection rates 
of X-ray, optical and radio afterglows. The mean flux 
densities and luminosities of radio afterglows were given 
in Frail (20051) . while Soderberg et al. (2006a)) compiled 
radio light curves t o compare GRB and sup e rnova l umi- 
nosities. Recently Ide Ugarte Postigo et al.l (|2011b[ ) has 
compiled the mm/submm data of GRBs and carried out 
the cumulative analysis in view of upcoming Atacama 
Large Millimeter Array (ALMA). Since the launch of 
the Swift satellite the number of GRBs with radio data 
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has doubled. Moreover, in the near future the contin- 
uum sensitivity of the Very Large Array^ (VLA; now 
Expanded Vary Large Array) - the primary telescope for 
radio afterglow follow-up - will increase by a factor of 5 
to 20 (depending on wavelength). It is therefore timely 
to pull together all the past radio data in one place and 
to use the catalog to define the average properties of ra- 
dio afterglows, search for trends in different sub-classes 
of GRBs (e.g. short vs long), search for differences in 
the Swift and pre-Swift sample, and to investigate if there 
are any correlations between the radio afterglows and the 
prompt or afterglow emission at shorter wavelengths. 

In [21 we present a summary of our radio sample and 
its properties. We also compare the detection statistics 
with respect to other wave bands. We discuss the de- 
tails and statistics of detection versus upper limits of our 
sample in ij3l In ij4]and [JH we investigate the distribu- 
tion and correlative properties of our radio sample with 
respect to the GRB prompt emission properties as well 
as the X-ray and optical afterglow properties. We plot 
the synthetic radio light curves and their dependence on 
various parameters in fJH We use these light curves to 
make some predictions for the Atacama Large Millimeter 
Array (ALMA) in [jTl while in [J8] we bring together to 
summarize the various results of the paper. 

We adopt a A-CDM cosmology throughout this paper 
with Hn = 71 km s" ^ Mpc-\ 17™ = 0.27 and = 0.73 
(jSpergel et al.l[2?)0l . 

2. THE RADIO SAMPLE AND STATISTICS 
2.1. Radio Sample 

Our compiled sample consists of 304 GRBs observed 
with the radio telescopes between January 1997 and Jan- 
uary 2011, along with the 2011 April 28 Fermi burst, 
GRB110428A. The sample consists of a total of 2995 
flux density measurements taken between the frequen- 
cies 0.6 GHz to 660 GHz bands spanning a time range 
from 0.026 to 1339 days. In Figure [U we plot the his- 
togram summarizing the distribution of the data in our 
catalog, as a function of frequency and time elapsed since 
the burst. A total of 1539 measurements were taken in 
8.5 GHz frequency band, while 657 measurements were 
taken in 5 GHz frequency band. 

The sample includes 33 Short-hard bursts (SHBs), 19 
X-ray flashes (XRFs) and 26 GRBs with possible su- 
pernova associations (SN/GRBs). The SN/GRB cat- 
egory includes all the GRBs with confirmed associa- 
tions as well as hkely supernova (SN) associations (see 
iHiorth fc Bloo^l20Tl . The SN/GRBs with known su- 
pernovae (SNe) associated with them in our sample are: 
GRB 980425 (SN 1998bw), GRB 011121 (SN 2001kc), 
GRB 021211 (SN 200211), GRB 030329 (SN 2003dh), 
GRB 031203 (SN 20031w), GRB 050525A (SN 2005nc), 
GRB 060218 (SN 2006aj), GRB 091127 (SN 2009nz) and 
GRB 101219B (SN 2010ma). The speclroscopically con- 
firmed SN associations among these are GRBs 980425, 
030329, 031203, 060218 and 091127. Of the 26 SN/GRBs 
associations, 17 and 9 were discovered in the pre-Swift 
and post-Swift epoch, respectively. Among the 33 SHBs, 

The Very Large Array is operated by the National Radio As- 
tronomy Observatory, a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, 
Inc. 
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Figure 1. Histograms summarizing the distribution of data of the 
radio selected sample catalog as a function of observing frequency 
(upper panel) and of the time elapsed since the burst (lower panel). 

28 were detected by the Swift alone. We also have one 
candidate galactic transient in our sample (GRB 070610, 
iKasliwal et all [20081 [Castro-Tirado et al.ll2008D . 

Most of the afterglows (270 in total) in our sample were 
observed as part of the VLA radio afterglow programs, 
whereas, 15 bursts were observed by the Expanded VLA 
(EVLA), and 19 southern bursts with the Australia Tele- 
scope Compact Array (ATCA). Among the total of 285 
VLA/EVLA bursts, 8 bursts were followed by both the 
VLA/EVLA and the ATCA; and 38 bursts were observed 
by the Westerbork Synthesis Radio Telescope (WSRT) in 
addition to the VLA/EVLA. A total of 19 VLA/EVLA 
bursts were followed by the Ryle telescope, and 11 bursts 
in our sample were observed by the Giant Metrewave Ra- 
dio Telescope (GMRT). The Very Long Baseline Array 
(VLBA) with its sub-milliarcsecond angular resolution 
was used to observe five bursts (see Table |6]). 

For the VLA observations, most of the afterglows were 
first observed in the 8.5 GHz frequency band, since it is 
the most sensitive VLA band. Once a GRB was detected 
in the 8.5 GHz band, we followed it at other VLA fre- 
quencies to measure the continuum spectrum. However, 
in cases of bright, nearby or high-redshift (high-z) GRBs 
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or GRBs with good amount of multiwaveband data, we 
carried out follow-up observations for few epochs even for 
the initially non-detected afterglows, depending upon its 
interest in other wavebands, VLA time availability fac- 
tors etc. The data were taken in standard intcrfcromct- 
ric mode for an integration time of typically 30 minutes 
(including the time on the calibrators and the GRB inte- 
gration time). The total bandwidth used was 100 MHz. 
The data was analysed using standard AIPS routines. 

In Table ini we list the comprehensive properties of our 
sample. In column 1, the GRB names are given. Col- 
umn 2 lists the missions under which these GRBs were 
discovered (please refer to the table notes for the abbre- 
viations of the mission names). We tabulate the J2000 
positions of the GRBs in columns 3 and 4. Columns 5, 
6 and 7 tabulate the detection statistics in X-ray, opti- 
cal and radio bands, respectively. Symbol 'Y' indicates 
a detection, while 'N' indicates a non-detection. The 
symbol 'X' means the burst was not observed in that 
particular band, whereas 'Y?' indicates that the obser- 
vation was made but a detection could not be confirmed. 
Column 8 tabulates all the radio telescopes used to ob- 
serve the corresponding burst. We report the Tgo du- 
rations (the time interval over which 90% of the total 
background subtracted counts are observed, with the in- 
terval s tarting at 5% of th e total counts have been ob- 
served (jKoshut et al.|[l995l )) in the observer's frame and 
the redshifts z in columns 9 and 10, respectively. Here 
Tgo durations are mentioned in the energy band of the 
specific detector which detected a GRB. For example, 
for Swift bursts, the Tgo durations are in 15-350 keV en- 
ergy range. Column 11 shows the 15-150 keV fluences 
(S'i5_i5o) of our sample and column 12 indicates the k- 
corrected isotropic bolometric 7-ray energy {E'^"!^ , energy 
range 1-10,000 keV in the rest frame). We report the X- 
ray fluxes at 11 hr {F^^'^) in 0.3-10 keV range in column 

13. Here the fluxes for the BeppoSAX bursts arc quoted 
in 1.6-10 keV range. We also report the R-band (0.7 fim) 
optical flux densities at 11 hr (F^^^, in fiJy) in column 

14. Column 15 tabulates the jet break times tj in cases 
where a clear jet break was seen and column 16 indicates 
the number density (n) of the circumburst medium. For 
the GRBs in which the number density of the medium 
was not known, we have assumed a density of 1 cm~'^ 
and have indicated them with "[1]". We estimate the 
coUimation angle dj (degrees) and the beaming-corrected 
bolometric energy (Et °i„) of GRBs with observed tj us- 
ing Frail ct al. (20QI|); iBloom et aH (|2003l ) and tabulate 
them in columns 17 and 18, respectively. Column 19 lists 
the references. Each cell has 8 references in the following 
order: Tgo, z, ^15-150, E^°l F|l^ F^l^ tj, and n. ^ If 
there is no data in a particular cell, the corresponding 
reference is indicated with []. However, a reference in- 
dicated by [] for a non-empty corresponding cell implies 
that that particular value has been calculated by us and 
this paper is the reference. 

2.2. Multiwaveband detection statistics of our sample 

In Table [U we summarize the detection statistics in ra- 
dio. X-ray and optical bands for all the bursts in our sam- 
ple. We subdivide the sample into pre-Swift and post- 
Swift bursts. This split was made to investigate possible 
selection biases. The p^e-Swift GRB event rate was low 




Figure 2. Venn diagrams illustrating the difTcrent relationships 
between the radio, optical, and X-ray afterglows for the entire sam- 
ple, the pic-Swift and the post-Swift sample. The symbols X, O 
and R indicates detections in X-ray, optical and radio bands, re- 
spectively. The combinations of these letters indicate detections 
in those respective bands. For example, RX indicates detection in 
radio and X-ray but not in optical, while ROX indicate afterglows 
seen in all three bands. 

enough tha t radio follow-up w as undertaken for nearly all 
events (see iFrail et all l2003| ) , while only a subset could 
be observed among the Swift bursts. 

Out of the 304 bursts, 123 bursts were observed in the 
j>TC-Swift epoch from 1997 until 2004. The remaining 
181 bursts were observed between 2005 and April 2011. 
A total of 95 bursts resulted in radio detections, while 
206 were non-detections. For three bursts an initial de- 
tection could not be confirmed. Thus the total detection 
rate for our sample in the radio band is 31%. The ra- 
dio detection statistics in pie-Swift and post-Swift sam- 
ples are 42/123 (34%) and 53/181 (29%), respectively. 
It is striking how similar the detection ratios are in both 
samples. In contrast, the X-ray detection rates increase 
from 42% to 93% after the launch of Swift. This is be- 
cause of the dedicated X-ray telescope (XRT) on-board 
Swift, which autonomously slewed to Burst Alert Tran- 
sient (BAT; onboard Swift) detected bursts. Note also 
that 51% of the bursts were not observed in the X-ray 
bands in the pre- Swift sample. The optical detection 
rates also increased from 48% to 75% between the pre- 
Swift and the post-Swift epoch. This too was due to 
the on-board UV/optical telescope (UVOT) which was 
autonomously slewed along with the XRT, and also due 
to the availability of rapid, well-localized positions for 
ground-based follow-up. We refer the reader to Table [T] 
for the detailed statistics of the multiwaveband observa- 
tions of our sample. 

In Figure [21 we plot a Venn diagram for our entire sam- 
ple illustrating the different relationships between the ra- 
dio, optical and X-ray observations of the afterglows. We 
also indicate this relationship for the pre- Swift and the 
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Table 1 

Multiwavcband statistics of the GRB afterglow sample 
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post-Swift sample separately. For the Venn diagrams, if 
a GRB did not have an observation in any band (i.e. 'X' 
in Tabic [Ij , we excluded them. We also excluded those 
events which had a unconfirmed detection in any bands 
(i.e. 'Y?'). Thus, we included a total of 226 bursts, 
out of which 10 bursts had no detection in any band. 
Among those 10 bursts with no afterglow detection, 6 arc 
SHBs (GRBs 021201, 020531, 050911, 051105A, 071112B 
and 090417A) and the remaining are normal long GRBs 
(GRBs 970111, 990217, 000615A and 071018). Hence, 
the final Venn diagram of the complete sample includes 
216 bursts. The X-ray detections are marked with X, op- 
tical detections are marked with O and radio detections 
are marked with R. ROX indicates the bursts which were 
detected in all three bands, whereas a combination of 
two letters indicates a joint detection in those respective 
wavelengths. For example, RX indicates bursts detected 
in radio and X-ray bands but not in optical. A single let- 
ter, X for example, indicates an afterglow detected in the 
X-ray band only. Our sample consists of 68 bursts which 
were detected in all three bands. There were 95 bursts 
which had detections in X-ray and optical but not in ra- 
dio. There was no burst which was detected in only the 
radio band (i.e. R), nor were there any bursts detected 
in only the radio and the optical bands without being 
detected in the X-ray band (i.e. RO). However, 6 bursts 
were detected in the radio and the X-ray bands but not 
in the optical. These form a subset of the dark bursts. 
Only one optical burst could not be detected in any other 
band. The number of only X-ray detected afterglows is 
46. In the post-Swift and the pre-Swift samples the X- 
ray only afterglows are 29 and 17, respectively. In the 
pre- Swift sample the GRBs detected in all three bands 
was only 18, rising to 50 in the post- Swift bursts. 

3. RADIO AFTERGLOW SAMPLE ANALYSIS 

3.1. Radio flux density distributions 

Since our data consists of detections as well as up- 
per limits, wc incorp orated the Kaplan-Meier Product 
Limit method (K-M; iFeigelson fc Nelsonlll985[ ) in deriv- 
ing the flux density distributions and mean flux density 
estimates of our sample. The K-M estimator was applied 
to a sample of bursts for which flux density measurements 
had been made at one frequency over a fixed time inter- 
val. For detections, we averaged multiple observations of 
the same GRB to obtain the mean flux density. This is 
perfectly acceptable since over the short periods the flux 
density variations at radio frequencies are dominated by 
interstellar scintillation. Table [5] summarizes our results 
for 5 and 8.5 GHz bands. Here, the mean estimates are 
given for only detections as well as by including the upper 
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Figure 3. The Kaplan-Meier probability distribution of the ra- 
dio flux densities at 8.5 GHz for a complete sample of 107 GRBs 
toward which measurements or upper limits were made between 5 
and 10 days. The distribution of the 51 detections is shown for 
comparison. 

hmits. The KM25, KM50 and KM75 arc the K-M method 
estimates of the minimum flux density above which 25%, 
50%, and 75 % of the total radio afterglows will lie. 

In Figure m we plot the Kaplan-Meier probability dis- 
tribution of the radio flux densities at 8.5 GHz for a com- 
plete sample of 107 GRBs toward which measurements 
or upper limits were made between 5 and 10 days. The 
distribution of the 51 detections is shown for compari- 
son. The overall detection rate in this time range is 48% 
{i.e., 51/107). The mean of the 51 detections is 275±34 
fxJy while the mean of the entire sample (upper limits 
included) is 170±20 /iJy. Half of all GRBs produce radio 
afterglows with flux densities at 8.5 GHz in excess of 92 

Before making any detailed comparisons, it is impor- 
tant to understand any possible biases that might ex- 
ist in these data. Specifically, we ask whether a radio 
non-detection is the result of an inadequate measurement 
(i.e., not observing at a correct time or with enough sen- 
sitivity) , or whether it is the result of an underlying phys- 
ical cause. In Figure |4l we plot the distributions for the 
radio 3-sigma upper limits as well as for the detections 
between 5-10 days. The upper limits peak in 100-150 
^Jy range, whereas the detections peak in 150-200 /iJy 
range with a tail extending to 1 mJy. As first noted by 
iFraill (I2005D , the difference in the upper limits and the de- 
tections is not highly significant (see also Table[3|). There 
is only about a factor of 50 difference between the radio 



A radio-selected sample of GRB afterglows 



5 



Table 2 

K-M estimates of the flux density distribution 
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219 
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203 
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302 ± 42 


185 


92 
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Note. — 

^ K-M estimates that 25% of all bursts 
^ K-M estimates that 50% of all bursts 
K-M estimates that 75% of all bursts 



will have radio afterglows above this value, 
will have radio afterglows above this value, 
will have radio afterglows above this value. 
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Figure 4. The filled histogram represents the distribution for the 
8.5 GHz 3-cr upper limits for GRIBs in our sample between days 
5-10. The hatched histogram shows the distribution for the detec- 
tions between days 5-10. The small range of peak flux densities 
between the radio detections and non-detections suggests that the 
detection fraction of radio afterglows is largely determined by the 
instrumental sensitivities. 

flux density of the brightest and the faintest detected 
cosmological bursts (i.e. z > 0.4). This low dynamic 
range stands in sharp contrast to the orders of magni- 
tude difference between the bright an d faint e vents seen 
in cither X-ray or afterglow samples (|Kann et al. 200^ 
iRacusin et aLlfeoilf ). This narrow flux density range sug- 
gests that the detection fraction of our radio afterglow 
sample is determined by the instrumental sensitivities. 

In Figure [Sj we plot all those GRBs which were never 
detected in the radio band. We also overlay one rare, 
bright GRB light curve GRB 980703 (red curve) and a 
more typical GRB 980329 (blue curve) . The plot shows 
that our upper limits do not cut off at a fixed flux density 
but are tightly clustered within one order of magnitude 
(between 30 //Jy to around 400 /xJy). As illustrated in 
Figure [SJ most of the observations would have been ca- 
pable of detecting a bright GRB 980703-like afterglow, 
while up to the first 10 days (when most of the afterglow 
observations take place), only around 50% observations 
would have been capable of detecting an average event 
like GRB 980329. This result reinforces our suggestion 
that radio afterglow searches are strongly sensitivity lim- 
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Figure 5. Upper limits at 8.5 GHz frequency band for all GRBs 
for which no afterglow was detected. The red line represents the 
light curve of a rare, bright event GRB 980703 and the blue line 
represents the light curve of a more typical event GRB 980329. The 
detection fraction of radio afterglows in the first 10 days certainly 
appears to be mainly limited by the sensitivity. Black dashed line 
indicates 3-sigma fiux density of the EVLA in its full capacity for 
a 30 min integration time. 

ited. In this figure, we also plot the 3-sigma flux density 
of EVLA for a 30min integration in its full capacity. 

3.2. k-corrected radio spectral luminosity 

In Figure [6l we plot the fc-corrected radio spectral 
luminosities with the 3-sigma upper limits. To con- 
vert the flux density into spectral luminosity L, we use 
L = ATTFd\/{\ + z), where F is the radio flux den- 
sity and dL is the luminosity distance corresponding to 
the redshift z. To incorporate ^-correction, we multiply 
this luminosity with a /c-corrcction factor of (1 -|- z)"'~^ 
(where a and jS are the time and frequency indices in 
F oc t"i/^). Thus the fc-corrected radio spectral luminos- 
ity IS L = AiTFd\{\-\-z)°'^'^^^ . Here we choose a = and 
/3 = 1/3 correspon ding to an optica lly thin, flat, post-jet- 
break hght curve (jFrail et al.ll2006l ). For the GRBs with 
unknown redshift, we assume the average redshifts of pre- 
Swift and Swift bursts, i.e. {z)pre = 1-3 for the pve-Swift 
burs ts and {z)post = 2.0 for the post-Swift bursts (see 
i )4.1[) . In our sample we have 147 GRBs with known red- 
shifts, while 157 GRBs don't have a measured redshift. 
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This figure, similar to Figures 2] and [SJ also shows that 
for most of the sample, the difference in the spectral lu- 
minosities between the detections and the upper limits 
is not very significant. The average luminosity for detec- 
tions is 1.1 X 10^^ erg s~^ Hz~^, whereas, the average lu- 
minosity 3-cr upper limit for non-detections is 6.4 x 10'^*' 
erg s~^ Hz~^, i.e., the difference between them is not 
even a factor of two. Figure |6] does show, however, the 
existence of the low-luminosity events which are several 
orders of magnitude fainter than the normal cosmolog- 
ical event s. W e will discuss these populations in more 
detail in g331 

3.3. Canonical light curves of radio afterglows 

In Figure [3 we plot light curves of all the detected 
radio afterglows at 8.5 GHz in their rest frames as well 
as the observer frame, excluding the GRB sub-classes 
(SHB, XRF and SNe/GRB) from which the low lumi- 
nosity events in Figure [5] originate. We average these 
data to produce a mean light curve for the long-duration, 
cosmological GRBs (LGRB) along with the 75% confi- 
dence bands. A canonical LGRB if observed from begin- 
ning to end will have a mean luminosity of ~ 2 x 10'^^ 
erg Hz~^ until about 3-6 d (10-20 d) in the rest- 
frame (observer-frame). After this time there is a grad- 
ual power- law decay with an index of approximately —1. 
The mean light curve shows two peaks (rest frame), one 
between day 0.1-0.2 and one around day 2, and a dip 
around day 1. The first peak is not significant since there 
are too few data points. The peak after day 2 is signif- 
icant but it is likely the same peak as identified in Fig- 
ure ini The dip near day 1 appears real and may signify 
a transition between two different emission components 
(reverse shock and forward shock). 

3.4. Bursts of different classes 

We explore the differences in the luminosity light 
curves of SHBs, XRFs and SN-GRBs. In Figure [51 
we plot the radio luminosity light curves for the differ- 
ent GRB classes in the 8.5 GHz band. The top plot 
is the light curve for SHBs. Our sample included 35 
SHBs but there are only two detections GRB 050724 and 
GRB 051221A. For 13 SHBs the redshift was not known. 
We averaged the redshifts of known-z SHBs and assigned 
this value ((z)=0.55) to the unknown-z SHBs. Here we 
also overlay the mean light curve of LGRBs obtained 
from the Figure [7l It is clear that SHBs are intrinsi- 
cally radio dim objects. They are more than an order 
of magnitude fainter than a LGRB. In the middle frame 
of Figure [5] we plot the light curves of the 19 XRFs in 
our sample which includes seven detections. Most XRFS 
are usually only about a factor of ten fainter than LGRB 
(red curve). The extreme outlier is GRB 060218, which 
is 4-5 orders of magnitude dimmer than a typical event. 
In the bottom frame of Figure [5] we selected only those 
9 GRBs for which the known SN associations is exist. 
The SNe/GRBs nearly fully populate this radio lumi- 
nosity plot. Some SNe/GRBs such as GRB 980425 and 
GRB 060218 are dim, while others like GRB 030329 are 
as bright 0jS cL typical LGRB. 

3.5. Peak of radio afterglow detections 

Out of the 95 radio-detected afterglows (see ^2.2^ . 63 
had radio light curves (i.e. 3 or more detections in a 



single radio band), whereas 32 bursts had less than 3 
detections. For the GRBs for which the light curves were 
available, we determined the peak flux density and the 
time of the peak at the VLA frequency bands (i.e. 1.4 
GHz, 4.9 GHz, 8.5 GHz, 15 GHz and 22.5 GHz bands) 
by fit ting the forward shock formula of the form (|Fraill 
l2005l) . 

IP f /—I f / 

This formula may not accurately represent the full 
complexity of the radio light curve evolution. However, 
it is good enough to determine the approximate values 
for the peak flux density Fm and the time of the peak 
tra- These fitted values of peak fiux density and peak 
times are provided in Table [31 All of the values quoted 
are in the observer's frame. In our sample there are 24 
radio afterglows for which we detect a small excess of 
radio emission at early times (t < 3 days) compared to 
the fit in Equation [l] We suspect that this excess may 
be the result of a separate reverse shock component. We 
will investigate this hypothesis in a separate paper. To 
avoid a possible reverse shock peak, we took the follow- 
ing precautions: If the light curve indicated two peaks, 
one before 3 days and one after that, we fit the later 
peak; If the best fit peak flux density was within flrst 
3 days, we exclude it as this could be the possible re- 
verse shock emission. The early-peaked light curve of 
GRB 060218 is included as the lone exception since it 
has a very low redshift, and hence the peak before day 
3 could still come from the forward shock. For the af- 
terglows which had less than 3 detections, we chose the 
value of highest fiux density. In the case of single detec- 
tion point at a particular frequency, we simply chose that 
value. There were 3 afterglows in 22.5 GHz band, 4 af- 
terglows in 15 GHz band, 7 afterglows each in 8.5 and 4.9 
GHz bands and 1 afterglow in 1.4 GHz band with single 
detection points. The second set of entries in Table[3|(i.e. 
after the horizontal line) represent the fiux density values 
taken directly from the data, and hence do not not have 
the best-fit errors. Rest frame peak times were obtained 
for the bursts with known redshift by using <„/(! + -z). 
We do not estimate the rest frame peak time for bursts 
with unknown redshifts. 



Table 3 

The radio peak flux densities and tlic time of peak for radio afterglows 



GRB 


Frcq. 


Peak Flux (F„) 


Peak time {tm ) 


t™/(l + ^) 




GHz 




days 


days 


970508 


1.43 


381 ± 19 


179.1 ± 7.6 


97.6 ± 4.1 


970508 


4.86 


780 ± 13 


57.6 ± 0.9 


31.4 ± 0.5 


970508 


8.46 


958 ± 11 


37.2 ± 0.4 


20.3 ± 0.2 


970828 


8.46 


144 ± 31 


7.8 ± 1.8 


4.0 ± 0.9 


980329 


4.86 


171 ± 14 


91.5 ± 11.3 


22.9 ± 2.8 


980329 


8.46 


332 ± 11 


33.5 ± 1.4 


8.4 ± 0.4 


980425 


1.38 


22260 ± 195 


47.1 ± 0.5 


46.7 ± 0.5 


980425 


2.5 


31621 ± 265 


32.7 ± 0.3 


32.4 ± 0.3 


980425 


4.8 


38362 ± 337 


18.3 ± 0.1 


18.1 ± 0.1 


980425 


8.64 


39360 ± 557 


12.7 ± 0.1 


12.6 ± 0.1 


980519 


4.86 


330 ± 47 


17.9 ± 2.4 




980519 


8.46 


205 ± 23 


12.6 ± 1.3 




980703 


1.43 


263 ± 81 


25.4 ± 5.5 


12.9 ± 2.8 


980703 


4.86 


1055 ± 30 


9.1 ± 0.2 


4.6 ± 0.1 


980703 


8.46 


1370 ± 30 


10 ± 0.2 


5.1 ± 0.1 


981226 


8.46 


137 ± 34 


8.2 ± 1.6 


3.9 ± 0.8 


990510 


8.46 


255 ± 34 


4.2 ± 0.5 


1.6 ± 0.2 
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Figure 6. The fc-corrected radio spectral luminosities at 8.5 GHz for radio-detected afterglows (black circles) versus non-detected 3-sigma 
luminosity upper limits (rod triangles) with respect to the rest frame time. The luminosity curve for the average cosmological burst varies 
over a small range, but there are also a number of low-luminosity events. 
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GRB 


Frcq. 


Peak Flux (F,„) 


Peak time {tm ) 


t^/{l + z) 


GRB 


Frcq. 


Peak Flux (F,„) 


Peak time (i^ ) 


t™/(l + 2) 




GHz 




days 


days 






GHz 




days 


days 


991208 


1.43 


263 ± 49 


8.9 ± 1.8 


5.2 ± 1 


1 


030329 


1.43 


2232 ± 30 


78.6 ± 2.1 


67.3 ± 1.8 


991208 


4.86 


845 ± 48 


12.8 ± 0.6 


7.5 ± 


4 


030329 


4.86 


10337 ± 33 


32.9 ± 0.1 


28.2 ± 0.1 


991208 


8.46 


1804 ± 24 


7.8 ± 0.1 


4.6 ± 


1 


030329 


8.46 


19567 ± 28 


17.3 ± 0.1 


14.8 ± 0.1 


000301C 


8.46 


520 ± 24 


14.1 ± 0.5 


4.6 ± 


2 


030329 


15 


34042 ± 84 


10.9 ± 0.1 


9.3 ± 0.1 


000418 


4.86 


897 ± 39 


27 ± 1 


12.7 ± 


5 


030329 


22.5 


49596 ± 66 


8.4 ± 0.1 


7.2 ± 0.1 


000418 


8.46 


1085 ± 22 


18.1 ± 0.3 


8.5 ± 


1 


030329 


43 


59318 ± 177 


5.8 ± 0.1 


5.0 ± 0.1 


000911 


8.46 


263 ± 33 


3.1 ± 0.3 


1.5 ± 


1 


030723 


8.46 


204 ± 40 


75 ± 10.6 




000926 


4.86 


460 ± 31 


16.9 ± 1.5 


5.6 ± 


5 


031203 


1.43 


929 ± 60 


65.5 ± 5.3 


59.3 ± 4.8 


000926 


8.46 


629 ± 24 


12.1 ± 0.5 


4.0 ± 


2 


031203 


4.86 


828 ± 28 


58.4 ± 2.1 


52.9 ± 1.9 


001018 


8.46 


590 ± 68 


4.7 ± 0.6 






031203 


8.46 


724 ± 19 


48 ± 1.3 


43.4 ± 1.2 


010222 


8.46 


93 ± 25 


16.8 ± 3.5 


6.8 ± 1 


4 


041219A 


4.86 


473 ± 28 


6.3 ± 0.9 




010921 


4.86 


161 ± 20 


31.5 ± 3.3 


21.7 ± 2 


3 


050416A 


4.86 


485 ± 36 


48.5 ± 3.5 


29.4 ± 2.1 


011030 


4.86 


121 ± 19 


19.4 ± 2.5 






050416A 


8.46 


373 ± 36 


49 ± 3.6 


29.7 ± 2.2 


011030 


8.46 


139 ± 9 


23.7 ± 1.4 






050509C 


8.46 


344 ± 45 


44.1 ± 5.7 




011121 


8.7 


655 ± 40 


8.1 ± 0.7 


5.9 ± 


5 


050603 


8.46 


377 ± 53 


14.1 ± 1.8 


3.7 ± 0.5 


011211 


8.46 


162 ± 13 


13.2 ± 1.6 


4.2 ± 


5 


050713B 


4.86 


424 ± 34 


21 ± 2.2 




020405 


8.46 


113 ± 17 


18.2 ± 3.7 


10.8 ± 2 


2 


050713B 


8.46 


343 ± 24 


14.1 ± 0.9 




020819B 


8.46 


291 ± 21 


12.2 ± 1.1 


8.7± 


8 


050820A 


8.46 


150 ± 31 


9.4 ± 1.9 


2.6 ± 0.5 


020903 


4.86 


782 ± 28 


36.7 ± 1.3 


29.4 ± 1 





050904 


8.46 


76 ± 14 


35.3 ± 1.5 


4.8 ± 0.2 


021004 


4.86 


470 ± 26 


32.2 ± 1.4 


9.7 ± 


4 


051022 


8.46 


268 ± 32 


5.2 ± 0.7 


2.9 ± 0.4 


021004 


8.46 


780 ± 23 


18.7 ± 0.5 


5.6 ± 


2 


060218 


4.86 


245 ± 50 


3.8 ± 0.7 


3.7 ± 0.7 


021004 


15 


1308 ± 260 


4.1 ± 0.8 


1.2 ± 


2 


060218 


8.46 


471 ± 83 


2 ± 0.3 


1.9 ± 0.3 


021004 


22.5 


1614 ± 52 


8.7 ± 0.3 


2.6 ± 


1 


070125 


8.46 


1028 ± 16 


84.1 ± 2 


33.0 ± 0.8 


021206 


4.86 


480 ± 69 


7 ±0.8 






070125 


15 


1250 ± 81 


18.6 ± 1.4 


7.3 ± 0.5 


030226 


8.46 


171 ± 23 


6.7 ± 1 


2.2 ± 


3 


070125 


22.5 


1778 ± 75 


13.6 ± 0.5 


5.3 ± 0.2 
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Figure 7. Upper panel: The radio light curves at 8.5 GHz for the 
long-duration, cosmological GRBs in the rest frame time. The red 
solid line represents the mean light curve. The pink shaded area 
represents the 75% confidence band. Lower panel shows the same 
plot for time in the observer's time frame. 
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GRB 


Frcq. 


Peak Flux (F,„) 


Peak time {tm ) 


t™/(l + z) 




GHz 




days 


days 


070612A 


4.86 


580 ± 20 


140.3 ± 5.8 


86.8 ± 3.6 


070612A 


8.46 


1028 ± 16 


84.1 ± 2 


52.0 ± 1.2 


071003 


8.46 


616 ± 57 


6.5 ± 0.5 


2.5 ± 0.2 


071010B 


4.86 


227 ± 114 


12.5 ± 4.4 


6.4 ± 2.3 


071010B 


8.46 


341 ± 41 


4.2 ± 0.5 


2.2 ± 0.3 


080603A 


8.46 


207 ± 26 


5.2 ± 1.8 


1.9 ± 0.7 


081221 


8.46 


174 ± 27 


8.9 ± 1.7 




090313 


8.46 


435 ± 22 


9.4 ±0.5 


2.1 ± 0.1 


090323 


8.46 


243 ± 13 


15.6 ± 1 


3.4 ± 0.2 


090328 


8.46 


686 ± 26 


16.1 ± 0.7 


9.3 ± 0.4 


090423 


8.46 


50 ± 11 


33.1 ± 6.3 


3.6 ± 0.7 


090424 


8.46 


236 ± 37 


5.2 ± 0.7 


3.4 ± 0.5 


090715B 


8.46 


191 ± 36 


9.2 ± 2.3 


2.3 ± 0.6 


090902B 


8.46 


84 ± 16 


14.1 ± 2.6 


4.9 ± 0.9 


091020 


8.46 


399 ± 21 


10.9 ± 0.6 


4.0 ± 0.2 


100414A 


8.46 


524 ± 19 


8 ± 0.3 


3.4 ± 0.1 


100418A 


4.86 


522 ± 83 


70.3 ± 8.1 


43.4 ± 5.0 


100418A 


8.46 


1218 ± 12 


47.6 ± 0.6 


29.4 ± 0.4 


100814A 


4.5 


496 ± 24 


13 ± 1 


5.3 ± 0.4 


100814A 


7.9 


613 ± 23 


10.4 ± 0.4 


4.3 ± 0.2 


990510 


4.86 


177 


9.2 


3.5 


991208 


15 


3100 


5.4 


3.2 
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Figure 8. A compilation of radio luminosity curves in the rest 
frame for three different populations of GRBs: short, hard bursts 
(top). X-ray hashes (middle), and supernova-GRB associations 
(bottom) . Here black filled circles represent the detections and the 
grey filled triangles represent the upper limits. The red solid line 
represents the mean light curve for the long-duration, cosmological 
sample from Figure [71 
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GRB 


Frcq. 


Peak Flux 


Peak time 










GHz 






days 




days 
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In Figure [9l we plot the histogram for rest frame peak 
times at 8.5 GHz band. It is clear that most of the ra- 
dio afterglows peak between 3 — 9 days in the rest frame. 
This could be due to the Vm passing the radio frequencies 
around this time. Another possibility is that jet breaks 
in GRBs occur mostly around this time (see tableland 
references therein) and radio stops rising after the occur- 
rence of jet break. In the Figure [HI we also plot the peak 
luminosity distribution of all the detected afterglows at 
all epochs in 8.5 GHz band. The luminosity distribution 
peaks in the range 10'^^ — 10^^ erg s~^ Hz~^. 



3.6. 1.^ GHz analysis 

While there is much less data at 1.4 GHz than at 8.5 
GHz, obtaining the afterglow properties at this frequency 
is still useful. There are a number of wide-field imag- 
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Figure 9. Upper panel: Histogram for tiic time-to-pcak distri- 
bution in the rest frame for radio afterglows at 8.5 GHz. Most of 
tlic afterglows peak between 3 — 6 days. Lower panel: Distribu- 
tion of peak radio luminosity at 8.5 GHz. Most of the bursts peak 
in the range 10'^^ — 10^^ erg s^^ Hz^^. 



ing instruments coming on line later this decade which 
are optimized to work at this frequency. This includes 
telescopes such as Australian Square Kilometre Array 
Pathfinder (ASKAP), the phased array feed Apertif on 
the Westerbork Synthesis Radio Telescope (WSRT) , and 
the South African MeerKAT array. 

A total of 310 measurements were made at this fre- 
quency for 55 GRBs. Out of 310 measurements, 125 re- 
sulted in detections corresponding to 12 GRBs. Thus in 
terms of number of GRBs detected, the detection rate is 
22%, while the detection rate in terms of measurements 
made is 40%. However, this detection rate is likely an 
overestimate since our 1.4 GHz sample is strongly biased. 
Observations at this frequency were not taken unless a 
bright radio afterglow had been detected at 8.5 GHz. 

In Figure 1101 we plot our non-detections as 3-f7 up- 
per limits with respect to the epoch since explosion. We 
also overlay 1.4 GHz hght curves of GRB 980703 (red 
curve) and GRB 970508 (blue curve). The plot shows 
that almost all the GRB upper limits are tightly clus- 
tered within one to two orders of magnitude (between 60 
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Figure 10. 3-cr upper limits at 1.4 GHz frequency band for all 
GRBs for which no afterglow was detected. The red line represent 
the light curve of event GRB 980703 and the blue line represents 
the light curve of event GRB 970508. Wc also plot the 1.4 GHz 
EVLA S-CT upper limit for a 30 minute integration time. 

^ Jy to around 3000 /i Jy) . As illustrated in Figure [TOl 
most of the observations would not have been capable 
of detecting afterglows from these representative GRBs. 
This figure suggests that radio searches are severely sen- 
sitivity limited at 1.4 GHz. With better sensitivity or 
more telescope time (see Figure [T0| . some improvement 
is expected detectable but as we show in ij6l the optimal 
search strategy is to use higher observing frequencies. 

In Figure [TTl we plot the /c-corrected 1.4 GHz radio 
spectral luminosities for the LGRBs and compare them 
with the 3-sigma upper limits for the non-detections. 
The figure again clearly shows that no difference in the 
spectral luminosities between the detections and the up- 
per limits, reinforcing the severe sensitivity limitation 
in resulting in various non-detections. The mean 1.4 
GHz luminosity for the detected bursts is 2 x 10'^° erg 
Hz~^, and the mean 3-cr upper limit for the non- 
detections is 10^^ erg s^^ Hz~^. Due to our strong selec- 
tion bias these numbers likely overestimate the luminos- 
ity of a typical GRB afterglow at this frequency. 

4. PARAMETER DISTRIBUTIONS 

4.1. Redshift distribution 

In our radio-selected sample of 304 GRBs we have 147 
GRBs with known redshifts. This includes two bursts for 
which only a redshift range is known (GRB 980326 with 
z = 0.9 - 1.1 and GRB 980329 with z = 2- 3.9. We have 
99 GRBs with redshifts available in the post- Swift epoch 
as compared to just 48 bursts with known redshifts in 
pre-Swift epoch. The remaining 157 bursts had no red- 
shift estimates, although 9 of them do have upper limits. 
In our sample, 59 bursts had z < 1 and 88 bursts had 
z > 1. GRB 090423 is the object wit h highest confirmed 
spectroscopic redshift of z = 8.26 (jTanvir et al.l 120091 : 
iSalvaterra et al.l [20091 ). whereas the highest photomet- 
ric redshift in our sain ple (z =9.4) is for GRB 090429B 
(jCucchiara et al.ll2011[ ) . The high est redshift in pre- Swif t 
sample is 4.5 for GRB 000131 () Andersen et al.l l2000t ). 
GRB 980425 is the object with the lowest redshift with 
z = 0.0085 (jTinnev et al.|[l998l) . In Figure [H we plot 
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Figure 11. The detected fc-corrected radio spectral luminosities 
at 1.4 GHz band (black circles) versus non-detected 3-sigma lumi- 
nosity upper limits (red triangles) with respect to the rest frame 
time since burst. The luminosity curve for the average cosmological 
burst is tightly clustered between 10'"' — lO^'^ erg s~^ Hz~^. 

the redshift distribution of our sample. We plot the pre- 
Swift and the post-Swift bursts for comparison. The fig- 
ures indicate that the post- Swift bursts are mostly found 
in the redshift range of 0.8 — 6.4 (71 bursts out of 99), 
whereas the maximum number of pre- Swift bursts lies in 
the redshift range of 0.4 — 3.2 (39 bursts out of 44). Here 
we also plot the redshift distribution of radio-detected 
and non-detected sample. Most of the non-detected ra- 
dio afterglows lie in range z =0.8-1.6 and the detected 
radio afterglows lie between redshift of 0.8-3.2. There 
is obviously no correlation between the redshift of the 
burst and its detectability in the radio band. Further 
support to the fact that the redshifts of the sample of 
radio detected GRBs and non-detected GRBs are not 
from different popul ations come from the Kolmogorov- 
bmirnov (K-S) test (H ress et al1ll992f ). Comparing the 
redshifts of the two sample gives P = 0.61, strengthening 
the fact that the two samples are not from the different 
populations. 

The average redshift of our sample is (z) total = 1.8. 
The average redshifts in pre- and post- Swift bursts are 
{z)pre = 1.3 and {z)post = 2.0, respectively. The mean 
redshift for the pre- Swift bursts is similar t o the value 
(z)=1.4 derived by iJakobsson et al.l (|2006l ). however, 
they measure (z)=2.8 for a post-Swift optically-selected 
sample. Instrumental sensitivity may explain the origin 
of this bias in our radio-selected sa mple to slightl y lower 
redshifts for post-Swift. However, iFvnbo et al.l ()2009l ) 
has measured the average redshift of Swift sample to be 
(z)=2.2, based on a sample of 146 Swift bursts, which 
is closer to mean redshift of our r adio selected sample. 
As first noted by iCiardi fc Loebl (|2000[ ) and shown by 
iFrail et ahl (|2006[ ). the detection rate of radio afterglows 
is largely insensitive to redshift. This is shown in Fig- 
ure 1131 where we plot peak radio flux densities at 8.5 
GHz band as a function of (1 -I- z). The red line shows 
the effect of the negative /c-correction factor or radio af- 
terglows, caused by the effects of spectral and temporal 
redshift, offsetting the diminution in distance. However, 
despite this effect we note that the mean peak flux den- 
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Figure 12. Upper panel: Redshift distributions for the post- 
Swift (light gray histogram) and the pic- Swift bursts (dark gray 
histogram). The average redshift of our sample is Zav = 1.8. The 
average redshifts in the pre- and post- Swift bursts are l^z^p^-^ — 1.3 
and (z)post = 2.0, respectively. The redshift distribution of the 
complete sample is shown with the hatched histogram. Lower 
panel: The redshift distribution of the detected (hatched his- 
togram) versus non-detected (filled histogram) radio afterglows for 
the complete sample. 

sity at z '^2.8 is close to the sensitivity limit of existing 
instruments. Figure 1131 however, shows that more sen- 
sitive telescopes, such as EVLA, will be able to detect 
higher-z bursts and enable us to test this hypothesis. In 
the mean time when discussing any distance-dependent 
effects for our sample we will adopt our derived redshifts 
rather than those in the literature. 

4.2. Fluence distribution 

Since the fluence of Swift bursts is available in the 
15-150 keV energy range, all the fiuences quoted in 
Table |6] are given in this range. We have converted 
fiuences from the non- Swift telescopes to this energy 
range (15-150 keV) ado pting the method explained in 
iNvsewander et all (|2009D . The fiuences quoted in Tabled 
are the fiuences in the observer's frame. In Figure we 
plot the fluence distribution of the entire sample as well 
as for the radio-detected and non-detected sample. Even 




1 10 

(Uz) 

Figure 13. Radio peak flux density versus (1 -|- z) plot for radio 
afterglows with known redshifts. The radio measurements are in 
the 8.5 GHz band (in observer's frame). Blue diamonds denote 
SNe/GRBs, while all the gray circles denote long-duration, cos- 
mological GRBs. Red and green lines are the symbolic lines (with 
random normalization) indicating the two different scenarios of the 
flux density scaling. The green dashed line indicates if the flux den- 
sity scales as simply the inverse square of the luminosity distance. 
The red thick line is the flux density scaling in the canonical after- 
glow model which includes a negative-fc correction effect, offsetting 
the diminution in distance. 3-o" EVLA sensitivity is also plotted to 
show the EVLA capability to go deeper in z. 

though both the samples peak in the range of 10^^ — 10^^ 
erg cm~^, 176 out of 206 (~85%) non-detected radio af- 
terglows have the fluence values < 10~^ erg cm~^, while 
86% of the radio detected, i.e. 82 out of 95 bursts, have 
fluence values > 10~^ erg cm~^. Thus there is a clear 
overall trend of non-detected radio afterglows favoring 
the lower values of fluence as compared to the higher flu- 
ence values favored by the detected. The K-S test also 
gives the P value of P = 2.61 x 10"'' supporting the fact 
that the radio-detected and non-radio-detected GRB flu- 
ences are being derived from the different populations. 

4.3. Energy distribution of the prompt GRB emission 

We calculated the isotropic-cquivalcnt 7-ray energies 
E^"l for all bursts with known redshifts. The isotropic 
energies given in Table [5] are /c— corrected bolometric en- 
ergies. We adopt the rest-frame bandpass in the energy 
range 1 keV-10 MeV. We quote all the original refer- 
ences from which the energy values are taken. However, 
for those bursts in which this energy was not known, we 
obtained the information about the best-fit model of the 
prompt GRB spectrum from the GCN circulars archive^. 
The prompt emis sion is best fit by either band function 
(|Band et al.lll993f ) or a cut- off power law f u nction (GPL) 
or a simple power law (see iRacusin et al.l ()2009D for de- 
tails). We use d these best fit para meters and adopted 
the method of IBloom et al.l (|2001D to estimate the k- 
corrected energy in 1 keV~10 MeV range. Since GRBs 
are coUimated events, the beaming-corrected energy is 
smaller than the isotropic-equivalent energy. To obtain 
this, we identified all the GRBs with known jet breaks 
{tj) from the literature. We also obtained all available 

* http://gcn.gsfc.nasa.gov/gcn3_archive.html 
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Figure 14. Upper Panel: Histogram of the fluence distribution 
for our entire sample. Lower panel: Fluence distribution of the 
radio-detected sample (hatched histogram) versus the non-detected 
sample (filled histogram). 
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Figure 15. The isotropic energy distribution of our complete sam- 
ple. Overlapped hatched histogram shows the beaming-corrected 
energy distribution. 



circumburst densities n from the literature. If a GRB 
had a known jet break but no available density measure- 
ment, we assumed density n = 1 cm~^ (we denote it as 
[1] in Table [6]). Once having values of t,- and rt, we est i- 
mated the coUimation angle 6j using (IFrail et al.ll200l . 
The beaming fraction = (1 — cos6'j) w was es- 
timated to calculate the beaming-corrected bolometric 



energies Ef°l^, where E\ 



bol 



bol 

ISO ■ 



In Figure [T5l we plot the isotropic energy ijf"^ dis- 
tribution available for 144 GRBs of our sample (46 
GRBs in pre- Swift epoch and 98 GRBs in post- Swift 
epoch). While maximum GRBs lie in the energy range 
of 10^^ — 10^'^ erg, the energy range spread of the sample 
is between 1.6 x 10''® — 4.1 x 10^^ erg, i.e. 7 orders of 
magnitude (this includes all the GRBs including SHBs, 
XRFs and SN-GRBs). This energy distribution and 
wide spread t o lower energies is characteristic of other 
large samples (|Kann "eFalllMoh . The fi gure also shows 
the beaming corrected energy distribution (hatched his- 
togram). There are only 48 GRBs in our sample (exclud- 
ing upper limits) which had beaming-corrected energies 
estimated. Even though most of the GRBs lie between 
energy range 10^° — 10^^ erg, the total spread in energy 
range is from 2.4 x 10*® - 1.4 x 10^^ erg. This distribu 



tion is similar to the previous studies fe.g. lRacusin et all 
[2009h . 

In Figure 1161 we plot the isotropic-equivalent and 
beaming corrected 7-ray energies for the radio-detected 
and non-detected sample. In case of isotropic ener- 
gies, the non-detected sample peaks at 10^^ — 10^^ erg, 
whereas, the detected sample peaks at 10^'^ — lO^'' erg 
range. While 60 out of 95 radio detected bursts have 
isotropic-equivalent > 10^"^ erg, only 9 out of 206 ra- 
dio non-detected bursts have energies above 10^'^ erg. 
For beaming corrected energies, the non-detected sam- 
ple peaks at 10"*^ — 10^° erg, whereas, the detected sam- 
ple peaks at 10^° — 10^^ erg range. This indicates that 
bolometric energies can indeed be a good indicator of 
the radio detectability of a burst. The K-S test gives 
the P values for the isotropic and beaming corrected en- 
ergies of the two samples to be P = 9.97 x 10"'' and 
P = 3.47 X 10"^, respectively (see Table 4). This im- 
plies that the two samples are inconsistent with being 
drawn from the same distribution and 7— ray energetics 
is indeed a good indicator of the radio detectabihty of a 
burst, with isotropic-equivalent energy being the stronger 
indicator. 

4.4. Distribution of X-ray and optical fluxes 

The X-ray fluxes {F^'^) quoted in Table [6] are in 
the energy range 0.3-10 keV. The X-ray fluxes of Bep- 
poSAX bursts are given in the energy range of 1.6-10 keV 
range. All X-ray fluxes are corrected for Galactic ex- 
tinction and we quote the X-ray fluxes at a time of 11 
hr after the burst. This is partly because most GRBs 
have measurements at this time but also because it is 
expected that the cooling frequency at this time will 
usually lie between the optical and X-ray frequencies 
(jNvsewander et al.ll2009f ). The afterglow flux wiU then 
be independent of the circumburst density for observing 
frequencies above the cooling frequency. We obtained 
most o f our X-ray fluxes at 11 hr from iGehrels et al.l 
(|2008f) : [de~Pasauale et al.l (|200l) : [Sakamoto et all (|2008l . 
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Figure 16. Upper panel shows the isotropic-equivalent 7-ray en- 
ergy distribution of the radio-detected sample (hatched histogram) 
versus the non-detected sample (gray filled histogram). Lower 
panel shows the same distribution for the beaming corrected ener- 
gies. 
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Figure 17. Upper panel: The distribution of the X-ray fluxes for 
our sample. Lower panel: This shows X-ray flux at 11 hr distribu- 
tion of the radio-detected sample (hatched histogram) versus the 
non-detected sample (gray filled histogram). 



|2011[ ). However, for the bursts in which the flux at 11 
hr was not available, we obtained the X-ra y light curve 
from the Swift XRT Lightcurve Repository (jEvans et al.l 
|2007| ) and derived the flux as close to 11 hr as possible 
(i.e. within ±0.5 hr). If the flux at ~ 11 hr was not 
available, we fit the X-ray light curve with a simple pow- 
erlaw and extrapolated the flux to 11 hr. This was done 
for 33 afterglows. Only in 5 cases, when the X-ray light 
curve was not available, did we extrapolated the avail- 
able X-ray flux to 11 hr using temporal index {ax, in 
Fx oc r^ ) ax = -1.17 for long b ursts and ax ^ -l.Tl 
for SHBs (jNvsewander et al.|[2009l ). 

Optical R band flux densities [F]^^) are also obtained 
at 11 hr in units of /xJy. For those bursts in which we 
could not find the desired values from the published refer- 
ences, wc used the GCN circulars and then extrapolated 
to the reported values at 11 hr using the best fit time de- 
cay index available. If a decay index was not available, 
we used a temporal slope of a_R = —0.85 (in Fr cx t"^) 
for long bursts and a^ = —0.68 for SHBs. The bursts for 
which an R-band magnitude was not available, we used 



Table 4 

P values for Kolmogorov-Smirnov test for the 
radio detected versus non-dctccted GRB sample 



Parameter P value 

Redshift 0.61 
Isotropic equivalent 7-ray energy 9.97 X 10^^ 

Beamed 7-ray energy 3.47 X 10""^ 

Fluence 2.61 X 10"^ 

X-ray flux (0.3-10 koV) 3.61 X lO"** 

Optical fl-band flux density 1.30 x 10^^ 



spectral index of /^j,; = —1 (in fj^ oc i/^" ) to convert the 
magnitudes into R band. Wc also corrected the flux den- 
sities for Galactic extinction using Sclilcgel et al. (1998). 
All the magnitudes were converted into /iJy using online 
NICMOS unit conversion tool^. 

In Figure [T71 we plot the X-ray flux distribution of our 
sample which is spread between 10~^^ — 10^^" erg cm~^ 
s^^. The majority of bursts have X-ray fluxes in the 

^ http:/ /www. stsci.edu/hst/nicmos/tools/conversion_form. html 
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Figure 18. Upper panel shows the distribution of the optical flux 
densities for our sample. Lower panel shows the optical flux den- 
sity distribution of the radio-detected sample (hatched histogram) 
versus the non-detected sample (gray flUed histogram). 
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Lower panel shows 
the X-ray flux at 11 hr of the radio-detected and non- 
detected samples. While the radio-detected sample peaks 
between 10~^^ — 10"^"'^ erg cm~^ range, the non- 
detected samples peaks at an order of magnitude less 
flux between the range 10""'^'^ — 10~^^ erg cm~^ s~^. In 
Figure [HI we plot the optical flux density distribution 
of our sample which ranges from 0.1 /iJy to 10** ^Jy. 
Most of the optical bursts lie in 10-100 fxJy range. In 
this figure, we also plot the optical R-band fiux density 
at 11 hr for the radio-detected and non-detected sample. 
The non-detected sample peaks at 1-10 /iJy, whereas, 
the detected sample peaks at 10-100 /iJy. 

On average, the X-ray fiuxes of bursts are an indica- 
tor of radio detectability. Likewise, the dust-extinction 
corrected optical magnitude of the burst is also a po- 
tential indicator of the radio detectability of the burst. 
The K-S test also supports this fact with P values of 
P = 3.61 X 10-6 and P = 1.30 x 10"^ for the X-ray and 
the optical fiuxes, respectively (Table ID 

5. CORRELATIVE PROPERTIES 
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Figure 19. Peak radio flux densities at 8.5 GHz band taken 
from Table \3\ on the X-axis and their corresponding fluences in 
15—150 kcV band are plotted on the Y-axis. Blue diamonds de- 
note SN/GRBs and green squares denote XRFs, while all the gray 
circles denote long cosmological GRBs. 

Here we proceed to look for correlations between the 
peak radio flux density (or luminosity) at 8.5 GHz band 
with various GRB properties (such as fluence, isotropic 
and bcaming-corrcctcd energy) and the X-ray and optical 
afterglow properties (redshift, flux density). We estimate 
the correlation in terms of the Pearson's correlation co- 
efficient, or R-index. The R-index varies between -1 to 
-|-1. A negative correlation coefficient indicates negative 
correlation between the quantities. An R-index close to 
means no correlation or a weak correlation. R-index val- 
ues of > 0.5 indicates a significant correlation between 
the parameters. We use only LGRBs to obtain the cor- 
relation coefficient to maintain homogeneity for this cor- 
relation analysis. 

In Figure [T9l we plot the 15-150 keV fiuences of radio 
detected GRBs against their corresponding peak radio 
flux densities at 8 GHz band from Table |3l There is no 
correlation between the two. The Pearson's correlation 
coefficient, R-index, between the two is 0.02. The figure 
also shows that while the spread in fluence is four or- 
ders of magnitude, radio peak fiux densities are tightly 
clustered. 

We also plot the isotropic 7-ray energies E'^"!^ against 
their respective peak radio luminosities In Figure [20l 
There in no correlation between the radio spectral lu- 
minosity and the isotropic energy of GRBs with R-index 
of 0.12.' 

In Figure [21] we plot the X-ray fluxes at 11 hr for 
radio-detected GRBs against their corresponding peak 
radio flux densities at 8.5 GHz band from Table [3l There 
is no correlation between the two with a R-index value 
of -0.05. We also plot the R-band optical flux densities 
at 11 hr for the radio-detected GRBs against their cor- 
responding peak radio flux densities at 8.5 GHz band 
from Table [3l The two parameters seem to be corre- 
lated. With a correlation index of 0.62, Fj^^'* is the only 
parameter correlated with the peak strength of a radio 
afterglow. Despite this correlation. Figure [5T] still illus- 
trates the clustering of the radio points over a narrow 
flux density range. 
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Figure 20. The peak radio luminosity at 8.5 GHz versus isotropic- 
equivalent 7-ray energies for various bursts. Here blue diamonds 
correspond to the SN/GRB events, green squares XRF events and 
gray circles are long cosmological GRBs. 



Table 5 

R-indcx correlation coefficient between various parameters 
of GRB afterglows 



Parameter 1 


Parameter 2 


R- index 


Peak radio flux density 


Fluonce 


0.02 


Peak radio luminosity 


Isotropic bol. energy 


0.12 


Peak radio flux density 


X-ray flux 


-0.05 


Peak radio flux density 


Optical flux 


0.62 



In Table [5l we tabulate the R-index correlation be- 
tween various GRB parameters. The correlation coeffi- 
cients vary between -1 to -|-1. The correlation coefficient 
values clearly indicate that the strength of the peak radio 
flux density does not depend upon the fluence, isotropic 
7-ray energy or the X-ray flux. Only the optical flux 
density at 11 hr has a positive correlation with the peak 
strength of the radio afterglows. 

6. SYNTHETIC RADIO LIGHT CURVES 

In order to compare the properties of our radio-sample 
with predictions from the basic afterglow theory, we plot 
synthetic light curves of GRB radio afterglows and de- 
termine their dependence on various afterglow param- 
eters such as density, kinetic energy, redshift and mi- 
croscopic shock parameters. In Figure we plot the 
radio a fterglow flux density a t the median GRB redshift 
z = 3 pakobsson et al.ll2006D and density n = 10 cm"'^ 
for different frequencies i.e. 1.4 GHz, 8.5 GHz, 35 GHz 
(EVLA bands) and 250 GHz (ALMA band). Here we 
fix other parameters such as the isotropic kinetic energy 
Ekeaso ~ 10^'^ erg, the beaming angle 9j = 0.2 rad, 
the electron energy density eg = 0.1, the magnetic en- 
ergy density cb = 1%, and the electron spectral index 
p = 2.2. The GRB parameters that we use here are 
;ood averages from the afterglow broadban d modeling 
Panaitescu fc Kumari[200Tt lYost et al.ll2003D . 

While not unique, the 8.5 GHz light curve in Figure [55] 
reproduces the average properties of our radio-selected 
sample derived in ^ The mean flux density (Table [2]), 
the time-to-peak and the overall timescale of the 
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Figure 21. Upper panel: Pealc radio flux densities at 8.5 GHz 
taken from Table |3] are on the X-axis and their corresponding X- 
ray fluxes at 11 hr in 0.3-10 keV band are plotted on the Y-axis. 
Blue diamonds denote SN/GR.Bs and green squares denote XRFs, 
while all the gray circles denote long cosmological GRBs. Lower 
panel: Peak radio flux density versus the R-band optical flux at 11 
hrs. 

decay is as seen in the data. The initial rise of the syn- 
thetic light curve at 8.5 GHz is steeper than the canonical 
long-duration event (Figure [7]), but as we have already 
noted (e.g. i j3.3p . this may be the result of an additional 
early emission component, possibly a reverse shock. The 
light curves in this figure are generated with only a for- 
ward shock component. 

With this close agreement between the model and the 
data, we can now investigate the effects of changing the 
observing frequency. As expected, there is a clear trend 
of increasing peak flux density and decreasing time-to- 
peak with increasing frequency. For illustration we com- 
pute the EVLA 3-ct sensitivities at 1.4 GHz, 8.5 GHz 
and 35 GHz bands and ALMA band-6 for a 1-hr inte- 
gration. These are 40 fiJy, 9 /iJy and 17/iJy, and 42/iJy, 
respectively. The challenges of detecting typical long- 
duration event events below 1.4 GHz are immediately 
apparent. For the above given set of parameters, which 
arc used to generate these light curves, EVLA, ASKAP, 
and WSRT/Apertif will not be able to detect radio af- 
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Figure 22. Plot of radio flux density light curves for a "stan- 
dard" GRB afterglow at various frequencies for a redshift z =3 
and density of 10 cm~^. See text for more details. 

terglows at 1.4 GHz unless they are observed for long 
integration times. Higher frequencies are clearly favored 
for detecting radio afterglows. We expect a large increase 
in the fraction of detectable radio afterglows with ALMA 
and the EVLA, should higher observing frequencies be 
used than in the past. 

In Figure [231 we derive the dependence of the radio 
afterglow light curve on various parameters, such as den- 
sity, microscopic parameters and kinetic energy. Since 
there is only a weak dependence on the electron energy 
density (cg) and electron spectral index (p), we do not 
show these plots here. 

We adopt the constant density model for the circum- 
burst medium to generate these synthetic light curves. 
The plots are made for a frequency of 8.5 GHz to enable 
comparison with our existing sample. As noted above, 
future EVLA observations would better to be conducted 
at higher frequencies. Increasing the centimeter radio 
frequency modifies the flux density and timescales some- 
what but it does not affect the dependencies that we see 
below. 

In the top left panel of Figure [53] we plot the synthetic 
light curves of radio afterglows for various density val- 
ues n at a redshift z = 3. We fix other parameters as 
ExEa^o = 10^3 erg, 6^ = 0.2 rad, = 0.1, eg = 1%, 
and p = 2.2. The radio afterglow brightness is a strong 
function of density. The flux density depends upon two 
competing effects. There is an initial increase in radio 
emission due to the enhanced synchrotron emission as 
the density increases but as the density continues to in- 
crease there is a reduction in the radio flux density due 
to the increasing synchrotron self-absorption effect. The 
figure shows that the radio afterglow is brightest for den- 
sities between n = 1 — 10 cm~'^; suggesting that existing 
radio samples are biased to a narrow range of circum- 
burst densities. At lower densities, the afterglow is in- 
trinsically weak, whereas at higher densities, synchrotron 
self-absorption effects are suppressing the radio afterglow 
strength for a long time. This may also explain why some 
of the bright GRBs are dim in radio band. They may 
need late-time observations in order to be detected. 

To determine the dependence of radio afterglow as a 



function of es, we plot the lightcurves for various val- 
ues of fractional magnetic energy density in the top right 
plot of Figure [23l We derive the light curves for a burst 
at a z = 3 and for circumburst density n = lOcm""^. 
We fix rest of the parameters as Ekeaso — 10'^^ erg, 
6j = 0.2 rad, = 0.1, and p = 2.2 as before. There is a 
clear trend. The higher the magnetic energy density, the 
brighter the radio afterglow peak. This is not surprising 
since the synchrotron emission is more efficient in high 
magnetic field. This trend is seen from n = 0.1 cm~'^ 
to n 10 cm~'^ but it does not continue to the higher 
density n = 100 cm""^. At higher densities the peak 
flux density does not increase dramatically with increas- 
ing eg, there is only a shift of the time-to-peak to later 
times. This is due to the effects of the large synchrotron 
absorption frequency which suppresses the emission. 

As expected, the radio afterglow strength strongly de- 
pends upon the kinetic energy of the burst (lower left 
panel of Figure I23p . The curve with an isotropic energy 
of 10^^ erg is identical to that in Figure (221 and we noted 
above that it provided a reasonable description of the av- 
erage properties of our sample. The beaming-corrected 
energy for this burst is 2 x 10^^ erg, a value that is at 
least consistent with the kinetic energy derived from ra- 
dio ca lorimetry (jBerger et al.ll2004bl : Ivan der Horst et al.l 
l2008h . 

In the lower right panel of the Figure [23l we determine 
the redshift dependence of radio afterglows for values of 
z =0.5. 1, 3, 5 and 8. Here we fix other parameters as 
Ekemo - 10^^ erg, 0j = 0.2 rad, = 0.1, eg = 1%, 
n = 10 cm"'^, and p — 2.2. Between the redshifts of 0.5 
to 3, the lightcurves get progressively fainter with de- 
creasing redshift, but beyond about z ^ 3, the negative-fc 
correction effects come into effect and the peak radio flux 
density falls only slightly with increasing redshift (see 
also HA.l^ . Thus, essentially all the bursts with kinetic 
enegies Ekeaso > 10^'^ erg should be easily detectable 
with less than 20-30 minutes exposure with the EVLA 
at any redshift. In Figure [25l we also plot some the the 
LGRBs with known redshifts between z = 0.8 to z = 8.3 
and show the capability of the EVLA to observe these 
bursts for a much longer duration. 

Finally in Figure [Mj we show flux density relationships 
for our standard afterglow model at X-ray, optical and 
radio wavelengths. These curves are not meant to be a 
complete exploration of the phase space but they rep- 
resent the broad behaviors to be expected with time. 
The correlations (or lack thereof) in ij5| can now be un- 
derstood. As we have already noted, our radio sample 
is sensitivity limited. Most 8.5 GHz measurements are 
tightly clustered between 30 ^Jy to around 400 /iJy with 
a mean of about 170 /iJy. In contrast, optical and X- 
ray measurements are sensitive to flux densities which 
vary over many more orders of magnitude. We can see 
from the figures that this limitation of the radio sample 
restricts the time range over which the correlations can 
be investigated. There is an additional complication in 
Figure [24| and that is the temporal evolution of the radio 
flux density is complex compared to the nearly mono- 
tonic behavior at X-ray and optical wavelengths. 

7. ALMA PREDICTIONS 

Our radio-selected catalog also contains millime- 
ter and submillimeter radio observations. There are 
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Figure 23. 8.5 GHz synthetic radio light curves for various values of density (top left), magnetic energy density (top right), isotropic 
kinetic energy (lowrer left) and redshift (lower right). EVLA sensitivity in this band for 1 hr integration time is 9 fi]y. 




Radio Flux Density at 8GHz (|ijy) Radio Fiux Density at SGHz (fjjy) R-band Opticai Flux Density (ijjy) 

Figure 24. Synthetic model radio flux density values at 8.5 GHz with respect to the corresponding R-band optical (left panel) and 1 keV 
X-ray (middle panel) flux densities. Time stamps are indicated from 0.1 to 10 00 days. For comparison, the right panel shows the plot of 
optical versus X-ray flux densities. The model is the same as that used in Fig. 1221 See text for more details. 
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Figure 25. 8.5 GHz radio light curve luminosity plots for some 
radio afterglows at different redshifts. We also plot the 1 hr EVLA 
3-0" sensitivity limit at a 2 = f , 3 and 8. 

186 mm/submm observations of 30 GRBs in our 
database so an analysis similar to the one we have 
done for the centimete r data is possible. However, 
Ide Ugarte Postigo et al.l (l2011bD have recently done a 
comprehensive analysis of all mm and submm data avail- 
able for GRB afterglows. It is more useful to simply 
summarize their findings and then compare their results 
with a typical GRB afterglow model (Figure [22|) and ex- 
plore the dependence of the radio millimeter afterglow 
brightness as a function of density, redshift, energy and 
es (e.g. m- 

The first mm/submm observations of GRBs were un- 
dertaken shortly after the discovery of the cm afterglow 
(iSmith et al.lll997at [Shepherd et aLlll998t iBremer et al.l 
1998[) . In subsequent compilations (jSmith et al.l 119991 
20051 ) it was well-established that the mm/submm detec- 



tion rate was strongly limited by sensitivity, with most 
detected events with fiux densities in excess of one or 
two milliJanskys. The situation has not improved sig- 
nificantly in subsequent years but the upcoming ALMA 
represents a jump of two orders of magnitude in contin- 
u um sensitivity. 

Ide Ugarte Postigo et al.l (|2011b[ ) report a detection 
rate of 25% for a sample of 102 mm/submm afterglows. 
Ide Ugarte Postigo et al.l ()2011bf ) use the detections and 
upper limits (assuming a Gaussian distribution) to make 
a rough estimate of the mean peak mm/submm fiux den- 
sity of 0.3 mJy. The average redshift of their sample is 
Zav = 1.99 and the average peak luminosity is io32 i±0-7 
erg s^^ Hz^^. 

Using typical GRB parameters, we showed in Fig- 
ure |22|) that a radio afterglow will be an order of mag- 
nitude brighter in the mm/submm (e.g. 250 GHz band) 
than in cm (e.g. 8.5 GHz). Support for this claim comes 
from average peak spectral luminosities which also differ 
by an order of magnitude ( §3.5p . 

We also can explore synthetic light curves identical to 
the cm curves from |j6]but applied to the ALMA band 6 
(211-275 GHz band) as a function of various parameters. 
In the upper left panel of Figure [211 we plot the synthetic 
radio light curves of GRB afterglows for various values of 
density n at a redshift z = 3. We fix other parameters as 



given in [JSj The mm/submm radio afterglow is a strong 
function of density. Since the effects of synchrotron self- 
absorption are weak in mm/submm bands, there is no 
competing effect to reduce the emission at higher densi- 
ties (unlike at cm bands). With its full collecting area, 
ALMA's 3-cr sensitivity for 1 hr integration time in band 
6 (211 -275 GHz) wiU be 42 fiJy (jde Ugarte Postigo et all 
I2011b[) . Thus ALMA should be able to detect all mm af- 
terglows for ri > 0.1 cm^'^. Not only is the mm/submm 
flux density brighter than the cm emission, it does not 
suffer from diffractive and refractive interstellar scintil- 
lation. 

To determine the dependence of radio afterglow bright- 
ness as a function of es, we plot the lightcurves for var- 
ious values of fractional magnetic energy density in the 
top right plot of Figure [551 We derive the light curves 
for a burst at a z = 3 and for circumburst density n ~ 10 
cm^"^. Here again we fix the rest of the parameters as 
above. We note that higher the magnetic energy density, 
the brighter the radio afterglow during its full evolution. 
Again this is in contrast to radio light curves at 8 GHz 
(iJB]), where higher was suppressing the optically thick 
flux density in high density medium. 

As expected, the mm afterglow strength also very 
strongly depends upon the kinetic energy of the burst 
(lower left panel of Figure [25]). ALMA will have some 
difficulties detecting a burst with Eke iso < 10^^. In 
the lower right panel of the Figure 1261 we determine 
the redshift dependence of radio afterglows for values 
of z =0.5, 1, 3, 5 and 8. Here we again fix other pa- 
rameters as Ekeaso = 10^"^ erg, 9j — 0.2 rad, ~ 0.1, 
eg = 1%, n = lOcm"'^, and p = 2.2. The trend is simi- 
lar to the cm but the mm/submm lacks the synchrotron 
self-absorption contribution to the negative-A: correction. 

8. DISCUSSION AND CONCLUSIONS 

In this paper we have formed a radio-selected sample 
of 304 GRBs based on 14 years of observations. In ad- 
dition to the radio observations, we have collected sup- 
plementary data including redshift, gamma-ray fluence, 
and optical and X-ray flux densities on all bursts. We 
have used these data to carry out the first comparative 
and correlative study of the prompt and afterglow emis- 
sion across the full electromagnetic spectrum. Previous 
studies (e.g . iGehrels et al.ll2008t iNvsewander et al.ll2009l: 



IKann et al.! 20101 ) have not included radio data in the 
analysis. 

We show that the fractional detection rate of radio 
afterglows is 31% ( tj2.2p . These detection statistics do 
not change substantially between the ■pie-Swift (42/123 
or 34%) and post-Swift (53/181 or 29%) samples. This 
is markedly different than the factor of two increase in 
the detection probability of optical and X-ray afterglows 
post-Swift. The benefits of Swift's on-board detection, 
autonomous slews, and the availability of rapid, well- 
localized positions are not realized by the slower evolving 
radio afterglows. 

Next we investigated the flux density distribution of 
the radio afterglows for cosmological GRBs (' ij3.5p . The 
mean flux density in the 8.5 GHz band between days 
2-8 of the 65 detections is 287±26 /iJy. But if we in- 
clude the non-detections (i.e. total of 157 8.5 GHz flux 
density points between days 2-8), using a Kaplan-Meier 
estimator, the mean is 162±15 /iJy. We also showed 




Figure 26. ALMA 250 GHz (Band 6) synthetic light curves for various values of density (top left), magnetic energy density (top right), 
isotropic kinetic energy (lower left) and redshift (lower right). ALMA band 6 sensitivity in its full capability will be 42 fi3y in 1 hr 
integration time. 



that range between the peak flux density of the brightest 
afterglow and the faintest afterglow spans only a nar- 
row range (a factor of 50), in sharp contrast to optical 
and X-ray afterglows which range over several orders of 
magnitude. Together this suggests that the radio radio 
afterglow searches are sensitivity limited and that this 
bias will be present in correlative studies (see below). 

Using the measured redshifts for our sample, we calcu- 
lated the mean spectral luminosity and the mean time- 
to-peak to be ^ 10'^^ erg s~^ Hz~^ and 3-6 days, re- 
spectively ( §3.5p . We investigated the luminosity dis- 
tribution of several sub-classes of GRBs including the 
long-duration, cosmological GRBs (LGRB), short-hard 
bursts (SHB), X-ray flashes (XRF) and GRBs with firm 
supernova associations (SNe/GRB). There arc a wide 
spread of luminosities for these different subclasses with 
LGRBs the brightest on average, while SHBs are nearly 
two orders of magnitude fainter than LGRBs, and XRFs 
and SNe/GRBs are ten times fainter than LGRBs, al- 
beit with large scatter ( ^3.4p . The best sampled radio 
light curves are those of the LGRB class with multi- 
ple data points that range from nearly 0.01 d to 1000 
d. The canonical LGRB will reach a peak luminosity 



of ^ 2 X lO'^^ erg s~^ Hz~^ about 3 to 6 days in the 
rest frame, and after about 10 to 20 days will undergo 
a power law decline with an index of order unity ( §3.3|) . 
There is some (not strong) evidence of two equally lu- 
minous emission components with a transition occurring 
at 1 d in the rest frame. The early component may be 
the result of the reverse shock as has been clai med previ- 
ously for several b ursts (.Kulkarni et al. 1999i: iFox et al.l 
[2003l : lBerger et al.ll2003briGhandra et al.ll2010[ ) while the 
later component is the standard forward shock which 
has been studied throiigh ex tensive broad-band modeling 
(jPanaitescu fc Kumaill2001[ ). 

Before any correlative analysis was carried out, we in- 
vestigated the high energy and non-radio afterglow prop- 
erties for our sample. We find that the distribution 
of gamma-ray fluence, duration, energy, and optical/X- 
ray flux densities for the 304 GRBs in our sample are 
not dissimilar those from previous studies. There is 
a difference in the redshift distribution ( ij4.ip . While 
the mean redshift for our pre-Swift sample agrees with 
past work, for our post-Swift sample we derive a lower 
value for the mean redshift (2.0) based on the radio se- 
lected sample compared to published values of z = 2.2 
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(iFvnbo et al.ll2009l] and z = 2.8 (jJakobsson et al.ll2006D . 
with iFvnbo et all ()2009D value closer to our mean Ow- 
ing to a negative fc-correction (jCiardi fc Loebll2000[ ). ra- 
dio afterglows are rather insensitive to redshift but we 
show that the mean radio flux density of high redshift 
GRBs are close to the sensitivity limit of a typical ob- 
servation made with existing instruments. This bias will 
have an effect on any distance-dependent correlations. 

We examined the observed GRB and afterglow param- 
eters for the radio sample (0 and we found that radio- 
detected GRBs had higher fluenccs, larger energies, and 
brighter X-ray and optical fluxes on average, compared to 
radio non-detections. There is evidently a statistical re- 
lationship between the detectability of a radio afterglow 
and these properties. Further correlations were explored 
between the peak radio flux density and gamma-ray flu- 
ence, optical flux density and X-ray flux density (|J5]). 
The only significant correlation that was found was be- 
tween the peak radio flux density and the optical flux 
density at 11 hr. 

The synthetic light curves Sj6] show that EVLA Ka 
band (26-40 GHz) and ALMA mm/submm bands are the 
most favorable bands for radio afterglow studies. We also 
show that the centimeter radio afterglow emission will be 
brightest for circumburst densities between n = 1 — 10 
cm~^. Outside of this range, the flux density will be weak 
either due to a low intrinsic emission strength (at lower 
densities) or due to increased synchrotron self- absorption 
(at higher densities). Within this density range there is 
a simple relationship of increasing strength of the cen- 
timeter radio emission and the fraction of the shock en- 
ergy in the magnetic field (es), but at higher densities 
synchrotron absorption suppresses the emission. The 
synthetic light curves also confirm what is already well- 
established and that is that the radio afterglow fiux den- 
sity is only weakly dependent on redshift at z > 2.5. 

We also make predictions at mm/submm wavelengths 
O: showing that ALMA should be able to detect 
most afterglows assuming typical parameters. The syn- 
chrotron self absorption effect docs not play at a major 
role in the mm/submm bands, thus ALMA should be 
able to detect afterglows very early on. ALMA should 
be able to provide a true distribution of the mm/submm 
afterglows, which is not well known. With an order of 
magnitude of continuum sensitivity of the EVLA and two 
orders of magnitude sensitivity of the ALMA, EVLA and 
ALMA will dominate GRB afterglow studies in the fu- 
ture. These telescopes will increase the temporal phase 
space and the redshift space over which GRBs can be 
observed. We illustrate this in Figure US] with the EVLA 
3-cr spectral luminosity sensitivity at redshifts of 1, 3 and 
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Xb 




Rd 


Radio 


^90 




z 


'^15-150° 






1 1 /i h 
R 


tj 








RofoponcoB'' 








.12000 


,12000 








Tel 








erg cm 


— A 


erg 


CGS 




d 




dcg. 


erg 




050813° 


s 




16:07 


57.00 


4-11:14:52.0 


Y 


N 


N 


VLA 


0.5 





720 


4 


40E 


08 


1.50E + 50 


0.0 












45, 32, 21, 8, 37, [],[],[] 


050814 


s 




17:36 


45.39 


+46:20:21.6 


Y 


Y 


N 


VLA 


151 


5 


300 


2 


OlE 


06 


6.00E + 52 


13.8 


2.1 


1.03 


[11 


3.05 


S.52E + 49 


45,37,21,8,37,37,42, [] 


050815 


s 




19:34 


23.15 


+09:08:47.5 


Y 


Y 


N 


VLA 


3 






7 


90E 


08 




0.9 












45, [], 45, [], 37, [],[],[] 


050820A 


s, 


T 


22:29 


38.11 


+ 19:33:37.1 


Y 


Y 


Y 


VLA 


240 


2 


615 


3 


44E 


06 


2.00E + 53 


221.2 


87.2 


7.35 


0.1 


5.07 


7.82E + 50 


8,33,21,8,37,21,42,10 


050824^ 


s 




00:48 


56.05 


+ 22:36:28.5 


Y 


Y 


Y 


VLA 


23 





830 


2 


66E 




1.50E + 51 


13.7 


9.9 










45,33,21,8,37,37,[],[] 


050904 


s 




00:54 


50.79 


+ 14:05:09.4 


Y 


Y 


Y 


VLA 


174 


6 


290 


4 


83E 


06 


1.30E + 54 


0.5 


20.0 




680 






45, 33, 45, 8, 37, 37,0,18 


050906° 


s 




03:31 


11.75 


-14:37:18.1 


X 


N 


N 


VLA 


0.3 






6 


OOE 


09 
















45, [], 45, [],[], [],[],[] 


050911° 


s 




00:54 


37.70 


-38:50:57.7 


N 


N 


N 


VLA 


16.2 





165 


3 


17E 


07 


2.69E+49 














45, 6, 45, 32, [],[],[],[] 


050915A 


s 




05:26 


44.80 


-28:00:59.3 


Y 


Y 


N 


VLA 


52 






8 


50E 


07 




5.1 


1.0 










45,[],21,[],37,37,[],[] 


050922B 


s 




00:23 


13.20 


-05:36:16.4 


Y 


N 


N 


VLA 


151 






2 


23E 


06 




20.0 












45,[1,45,[1,37,[1,[],[1 


050922C 


H 


s 


21:09 


33.30 


-08:45:27.5 


Y 


Y 


Y 


VLA 


5 


2 


199 


1 


62E 


06 


3.90E + 52 


17.2 


14.3 










45,33,21,S,37,21,[],[) 


051006 


s 




07:23 


13.52 


+09:30:24.5 


Y 


Y 


N 


VLA 


35 






1 


34E 


06 




0.6 


11.4 










45,[1,45,[1,37,37,[],[1 


051008 


s, 


T 


13:31 


29.30 


+42:05:59.0 


Y 


Y 


N 


VLA 


280 






5 


09E 


06 




20.6 


0.4 










36,[1,21,[],37,37,[],[1 


051016A 


s 




08: 1 1 


16.30 


-18:17:49.2 


Y 


N 


N 


VLA 


23 






8 


37E 


07 




3.1 












45,[1,45,[1,37,[1,[],[1 


051016B 


s 




08:48 


27.60 


+ 13:39:25.5 


Y 


Y 


N 


VLA 


4 





936 


^ 


70E 


07 


3.70E + 50 


15.5 


4.9 


1.56 


[1] 


10.49 


6.20E + 48 


36, 33, 21, 8, 37, 37, 42, [] 


051021A° 


H 




01:56 


36.39 


+09:04:03.7 


Y 


Y 


N 


VLA 


27 






3 


78E 


06 




13.7 


6.2 










26,[1,(],[1,37,37,[],[] 


051022 


H 


T 


23:56 


04.10 


+ 19:36:24.2 


Y 


Y 


Y 


VLA'^ 


200 





809 


8 


25E 


05 


6.30E + 53 


426.5 


<0.4 










26,9,[],3,37,37,[],[] 


051103° 


T 




09:52 


34.08 


+ 68:50:42.0 


X 


N 


N 


VLA 


0.2 






8 


66E 


06 
















26, [],[], [],[], [],[],[] 


051105A° 


S 




17:41 


03.28 


+ 34:59:03.6 


N 


N 


N 


VLA 


0.09 






2 


20E 


08 
















45, [1,45, [],[], [],[],[] 


051109A 


s, 


T 


22:01 


15.31 


+40:49:23.3 


Y 


Y 


Y 


VLA* 


37 


2 


346 


2 


20E 


06 


2.30E + 52 


99.7 


39.6 


0.92 


[1] 


4.18 


6.13E+49 


45, 33, 21, 8, 37, 21, 42, [] 


051109BP 


s 




23:01 


50.21 


+ 38:40:46.0 


Y 


N 


N 


VLA 


14 





080 


2 


56E 


07 


3.60E+48 


3.9 












45, 37, 45, 8, 37, [],[],[] 


051111 


s 




23:12 


33.36 


+ 18:22:29.5 


Y 


Y 


Y 


VLA 


46 


1 


550 


4 


08E 


06 


6.00E + 52 


3.4 


41.6 










45,33,45,8,37,37,[],[) 


051114° 


s 




15:05 


03.84 


+ 60:09:21.6 


N 


X 


N 


VLA 








1 


28E 


















(1, [1,46, [],[], [],[],[] 


051117A 


3 




15: 13 


34.09 


+ 30:52:12.7 


Y 


Y 


N 


VLA 


136 






4 


34E 


06 




11.3 


6.9 










45,[],21,[],37,37,[],[] 


051117B 


s 




05:40 


43.00 


-19:16:26.5 


Y 


Y 


N 


VLA 


9 






1 


77E 


07 




0.2 


<12.4 










45,[],45,[],37,37,[],[] 


051211B 


I 




23:02 


41.57 


+ 55:04:51.5 


Y 


Y 


Y 


VLA 


80 






1 


90E 


06 




2.9 


135.8 










55, [],[],[], 59, 37,[],[] 


051221A° 




T 


21:54 


48.63 


+ 16:53:27.2 


Y 


Y 


Y 


VLA'^ 


1.4 





547 


1 


15E 


06 


2.80E + 51 


9.9 


7.5 


4.07 


0.001 


5.35 


1.22E+49 


32,32,21,8,37,37,42,49 


051221B 


s' 




20:49 


35.10 


+ 53:02:12.2 


Y 


Y 


N 


VLA 


40 






9 


17E 


07 






<21.4 










45,[],45,[],37,37,[],[] 


051227° 


s 




08:20 


58.11 


+ 31:55:31.9 


Y 


Y 


N 


VLA 


115 





714 


6 


99E 


07 


8.00E + 50 


0.9 


0.5 










45,26,21,8,37,37,[],[] 


060105 


s, 


T 


19:50 


00.60 


+46:20:58.3 


Y 


N 


N 


VLA 


54 






1 


76E 


05 




37.8 












45, [1,45, [1,37, [],[],[] 


060108 


s 




09:48 


01.98 


+ 31:55:08.6 


Y 


Y 


N 


VLA 


14 


<2.8 


3 


69E 


07 




9.8 


3.0 










36,9,21,[],37,21,[],[] 


060110 


s 




04:50 


57.85 


+ 28:25:55.7 


Y 


Y 


N 


VLA 


26 






1 


57E 


06 




21.2 


59.4 










45,[1,45,[1,37,37,[],[1 


060116 


s 




05:38 


46.28 


-05:26:13.1 


Y 


Y 


Y 


VLA 


106 






2 


41E 


06 




5.8 


6.9 










45,[1,45,[1,37,37,[],[1 


060123 


s 




11:58 


47.66 


+45:30:50.6 


Y 


N 


N 


VLA 


900 


1 


099 


3 


OOE 


07 




27.6 












26, 26, 37, [1,37, [],[],[] 


060124 


H 


s 


05:08 


25.50 


+ 69:44:26.0 


Y 


Y 


N 


VLA 


800 


2 


297 


4 


OlE 


07 


7.00E + 51 


287.3 


68.8 


0.68 


[1] 


4.36 


2.02E+49 


36, 33, 21, 8, 37, 21, 42, [] 


060206 


s 




13:31 


43.44 


+ 35:03:03.2 


Y 


Y 


N 


VLA'^ 


8 


4 


050 


8 


31E 


07 


4.07E + 52 


20.8 


82.7 










45,33,21,33,37,21,[],[] 


060210 


s 




03:50 


57.37 


+ 27:01:34.4 


Y 


Y 


N 


VLA 


255 


3 


910 


7 


66E 


06 


4.20E + 53 


85.6 


1.3 










36,33,21,8,37,37,[],[1 


060213 


T 




09:26 


24. 14 


-09:06:57.4 


Y 


Y? 


N 


VLA 


60 






1 


60E 


04 
















26,[1,37,[1,[1,[1,[],[1 


060218^'P 


s 




03:21 


39.68 


+ 16:52:01.8 


Y 


Y 


Y 


VLA*'^ 


128 





033 


1 


57E 


06 


2.90E+48 


48.8 


1082.2 


>22 


5 


>80.26 


>2.85E+48 


8,33,45,8,37,37,50,50 


060313° 


s, 


T 


04:26 


28.40 


-10:50:40.1 


Y 


Y 


N 


VLA 


0.7 






1 


13E 


06 




6.6 


9.9 










32,[],21,32,37,37,[],[1 


060418 


s. 


T 


15:45 


42.40 


-03:38:22.8 


Y 


Y 


Y 


VLA 


103 


1 


490 


8 


33E 


06 


l.OOE + 53 


9.4 


25.0 










45,33,21,S,37,21,[],[] 


060421 


s 




22:54 


32.63 


+ 62:43:50.1 


Y 


N 


N 


VLA 


12 






1 


25E 






3.9 












45,[1,45,[1,37,[1,[],[1 


060428BP 


s 




15:41 


25.63 


+62:01:30.3 


Y 


Y 


N 


VLA 


58 






8 


23E 


07 




3.4 


6.3 










45,[],21,[],37,37,[],[1 


060502B° 


s 




18:35 


45.89 


+ 52:37:56.2 


Y 


Y 


N 


VLA 


0.13 





287 


4 


OOE 


08 


3.00E+49 


0.1 


<1 










45,32,21,S,37,37,[],[1 


060505° 


s 




22:07 


04.50 


-27:49:57.8 


Y 


Y 


N 


VLA 


4 





089 


9 


44E 


07 


4.37E+49 


8.1 


2.2 










32, 32, 45, 32, 37, 37, [],[1 


060522 


s 




21:31 


44.80 


+02:53:10.4 


Y 


Y 


N 


VLA 


71 


5 


110 


1 


14E 


06 


7.00E + 52 


3.2 


3.7 










45,36,45,S,37,37,[],[1 


060604 


s 




22:28 


55.01 


-10:54:55.8 


Y 


Y 


N 


VLA 


95 


2 


680 


4 


02E 


07 


4.37E + 51 


16.4 


17.2 










45,9,21,36,37,21,[],[] 


060605 


s 




21:28 


37.32 


-06:03:31.3 


Y 


Y 


N 


VLA 


79 


3 


773 


6 


97E 


07 


2.50E + 52 


6.1 


17.8 


0.17 


[11 


1.92 


1.41E+49 


45, 33, 21, 8, 37, 21, 42, [] 


060707 


s 




23:48 


19.00 


-17:54:17.0 


Y 


Y 


N 


VLA 


66 


3 


430 


1 


OOE 


06 


6.10E + 52 


25.0 


8.1 


12.26 


[1] 


8.80 


7.20E+50 


45, 33, 21, 8, 37, 37, 42, [] 


060719 


s 




01:13 


43.57 


-48:22:55.0 


Y 


Y 


N 


ATCA 


67 


<4.6 


1 


50E 


06 




21.4 


<0.1 










45,20,21,[],37,37,[],[] 


060801° 


s 




14:12 


01.31 


+ 16:58:54.0 


Y 


N 


N 


VLA^ 


0.5 


1 


131 


8 


OOE 


OS 


3.09E + 51 


0.1 












32,32,21,32,37,[1,[],[] 


060825 


s 




01:12 


29.00 


+ 55:47:46.5 


Y 


N 


N 


VLA 


8 






9 


63E 


07 




3.6 












45,[1,45,[1,37,[1,(],[1 


060908 


s, 


T 


02:07 


18.36 


+00:20:31.2 


Y 


Y 


N 


VLA 


19 


1 


884 


2 


80E 


06 


7.00E + 52 


5.9 


0.8 










36,20,45,8,37,37,[],[1 


060912A 


s. 


T 


00:21 


08.16 


+ 20:58:17.8 


Y 


Y 


N 


VLA 







937 


8 


69E 


07 


S.OOE + 51 


4.9 


10.3 










36,38,45,S,37,37,[],[] 


060923A 


s 




16:58 


28.15 


+ 12:21:38.9 


Y 


Y 


N 


VLA 


52 






8 


69E 


07 




10.8 


0.3 










45,[1,21,[],37,37,[],[1 


060923C 


s 




23:04 


28.49 


+03:55:26.6 


Y 


Y 


N 


VLA 


76 






1 


58E 


06 




3.6 


0.8 










45,[1,21,[],37,37,[1,[] 


060926 


s 




17:35 


43.66 


+ 13:02:18.6 


Y 


Y 


N 


VLA 


8 


3 


209 


2 


19E 


07 


l.OOE + 52 


2.5 


3.3 










45,20,45,S,37,37,[],[1 


061028 


s 




06:28 


54.66 


+46:17:57.0 


Y 


N 


N 


VLA 


106 





760 




66E 


07 


2.29E + 51 


0.8 












45,33,45,33,37,[1,[],[1 


061121 


s, 


T 


09:48 


54.58 


-13:11:42.7 


Y 


Y 


Y 


VLA'^'^ 


81 


1 


315 




37E 


05 


1.90E + 53 


115.1 


27.4 










36,33,21,8,37,21,[],[] 


061126 


s, 


T 


05:46 


24.43 


+ 64:12:38.8 


Y 


Y 


N 


VLA 


71 


1 


159 


6 


77E 


06 


S.OOE + 52 


44.4 


8.3 










45,33,21,8,37,21,[],[] 


061210° 


s 




09:38 


05.04 


+ 15:36:46.8 


Y 


N 


N 


VLA 


85 





410 


1 


HE 


06 


9.00E + 50 














45,32,45,S,[],[1,[],[1 


061222A 


s, 


T 


23:53 


03.42 


+46:31:58.6 


Y 


Y 


Y 


VLA 


72 


2 


088 


7 


99E 


06 


1.03E + 53 


282.1 


0.3 


0.77 


[1] 


2.39 


1.31E+51 


36,57,21,[],37,37,42,[] 


061222B 


s 




07:01 


24.60 


-25:51:36.0 


Y 


Y 


N 


VLA 


40 


3 


355 


2 


24E 


06 


S.OOE + 52 


3.1 


1.2 










45,37,45,8,37,37,[],[1 


070125 


T 








+ 31:09:04.1 


Y 


Y 


Y 


VLA^'^ 


60 


1 


548 






05 


9.55E + 53 


37.8 


207.7 


3.8 


42 


7.90 


9.07E+51 


26, 33, [1,33,37,37, 11, 11 


070220 


s 




02:19 


06.83 


+68:48:16.1 


Y 


N 


N 


VLA 


129 






1 


04E 


05 




10.6 












45,[1,21,[],37,[1,[1,[1 


070223 


s 




10:13 


48.39 


+43:08:00.7 


Y 


Y 


N 


VLA 


89 






1 


70E 


06 




6.9 


0.6 










45,[1,21,[],37,37,[],[1 


070306 


s 




09:52 


23.38 


+ 10:28:55.2 


Y 


Y 


N 


VLA 


210 


1 


497 


5 


38E 


06 


6.00E + 52 


16.1 


5.0 


1.33 


[11 


4.75 


2.07E+50 


45, 20, 21, 8, 37, 37, 42, [] 


070311 


I 




05:50 


08.42 


+03:22:30.0 


Y 


Y 


N 


VLA 


50 






1 


OOE 


06 




13.5 


8.1 










55,[1,[],[1,37,37,[],[] 


070429B° 


s 




21:52 


03.82 


-38:49:42.2 


Y 


Y 


N 


VLA 


0.5 





902 


6 


30E 


08 


1.35E + 50 


0.5 


<0.3 










32, 32, 45, 32, 37, 37, [],[1 


070518 


s 




16:56 


47.50 


+ 55:17:50.6 


Y 


Y 


N 


VLA 


5.5 






1 


62E 


07 




3.0 


3.3 










45,[1,21,[],37,37,[],[1 


0706101 


s 




19:55 


09.63 


+ 26:14:05.6 


Y 


Y 


N 


VLA^ 


10 






2 


40E 


07 






0.9 










26,[1,26,[1,[],26,[],[1 


070612A 


s 




08:05 


29.61 


+ 37:16:15.1 


X 


Y 


Y 


VLA^'*''^ 


369 





617 


1 


OOE 


05 


9.12E + 51 




329.8 










45,33,45,33,[],37,(],[] 


070714B° 


s, 


T 


03:51 


22.23 


+ 28:17:50.8 


Y 


Y 


N 


VLA 


3 





923 


7 


25E 


07 


l.lOE + 52 


1.4 


2.9 










32, 25, 46, 32, 37, 37, [],[1 


070721B 


s 




02:12 


32.97 


-02:11:40.4 


Y 


Y 


N 


VLA 


32 


3 


630 


3 


OOE 


06 


3.13E + 53 


5.2 


0.6 










26,20,21,S,37,37,[],[1 


070724A° 


s 




01:51 


13.96 


-18:35:40.1 


Y 


Y 


N 


VLA 


0.4 





457 


3 


OOE 


OS 


2.45E+49 


2.4 












32,32,21,32,37,[1,[],[1 


070724B 


sA 


> T 


01:09 


56.34 


+ 57:40:34.3 


Y 


N 


N 


VLA 


50 






1 


13E 


05 




84.8 












26,[1,[],[1,37,[],[],[] 


070729° 


s, 


T 


03:45 


15.32 


-39:19:26.4 


Y 


N 


N 


VLA 


0.9 






1 


OOE 


07 




0.2 












32,[],21,32,37,[1,[],[1 


070923° 


s 




12:18 


29.52 


-38:17:38.4 


X 


X 


N 


VLA 


0.2 






4 


92E 


OS 
















46,[1,46,[1,[],[1,[],[1 
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.12000 J2000 Tel s erg cm~^ erg CGS ^iJy 



071003 


s, 


T 


20 


07 


24 


22 


4-10:56:50.0 


Y 


Y 


Y 


VLA 


148 


1.604 


8 


30E 


06 


3.24E + 53 


55.7 


180.1 










46,33,37,33,37,37,[],[] 


071010A 


S 




19 


12 


14 


62 


-32:24:07.2 


Y 


Y 


N 


VLA 


6 


0.985 


2 


OOE 


07 


1.32E + 51 


34.9 


43.0 










46,33,37,33,37,37,[],[] 


071010B 


s, 


T 


10 


02 


09 


26 


+45:43:50.3 


Y 


Y 


Y 


VLA 


36 


0.947 


4 


40E 


06 


2.60E + 52 


45.6 


48.3 










46,37,37,46,37,37,[],[] 


071011 


S 




00 


33 


32 


73 


+ 61:08:04.6 


Y 


Y 


N 


VLA 


81 




2 


20E 


OS 




46.0 


26.6 


5.49 


[1] 






46,[],37,[],37,37,42,[] 


071018 


S 




10 


58 


44 


40 


+ 53:49:19.2 


N 


N 


N 


VLA 


288 




1 


OOE 


06 
















46, [], 37, [],[], [],[],[] 
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